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ABSTRACT 



■ Aims. The accretion of stars onto the central supermassive black hole at the center of the Milky Way is predicted to generate large 

fluxes of subrelativistic ions in the Galactic center region. We analyze the intensity, shape and spatial distribution of de-excitation 
gamma-ray lines produced by nuclear interactions of these energetic particles with the ambient medium. 

Methods. We first estimate the amount and mean kinetic energy of particles released from the central black hole during star disruption. 
O !■ We then calculate from a kinetic equation the energy and spatial distributions of these particles in the Galactic center region. These 

I particle distributions are then used to derive the characteristics of the main nuclear interaction gamma-ray lines. 

Results. Because the time period of star capture by the supermassive black hole is expected to be shorter than the lifetime of the ejected 
fast particles against Coulomb losses, the gamma-ray emission is predicted to be stationary. We find that the nuclear de-excitation 
lines should be emitted from a region of maximum 5° angular radius. The total gamma-ray line flux below 8 MeV is calculated to 
be «10~ 4 photons cirT 2 s -1 . The most promising lines for detection are those at 4.44 and ~6.2 MeV, with a predicted flux in each 
line of «10~ 5 photons cirr 2 s~'. Unfortunately, it is unlikely that this emission can be detected with the INTEGRAL observatory. But 
the predicted line intensities appear to be within reach of future gamma-ray space instruments. A future detection of de-excitation 
gamma-ray lines from the Galactic center region would provide unique information on the high-energy processes induced by the 
f"^ ' central supermassive black hole and the physical conditions of the emitting region. 
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1 ■ Introduction few months a circular transient accretion disk around the black 
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hole ( Ulmerl[l99 9i). The time-dependent accretion rate follows a 



3\ . Massive black holes (MBH) at galactic centers are sources of A k r ^ mjdt K r5/3 with an initial accretion higher 

^ ■ high energetic activity. X-ray observations of these sources re- than thg Eddington rate 

>vealed flares of hard X-ray photons releasing a maximum power . . . , , J c , . .. 

„ , n AA _i a iii A total tidal disruption of a star occurs when the penetration 

of 10 ergs . These flares are supposed to be due to processes ,_i , . , • tU c j- 

. b • r i . • i. . parameter b » 1, where b is the ratio of r„ - the penapse dis- 

of accretion and tidal disruption of stars by the massive black f , ,. . f , . us . K_, , J. f, r™. 

»j ■ , , . rr: riT^^^^ r , . . tance (distance of closest approach) to the tidal radius Rj. The 

holes (see, e.g., Sirota et al., 2005, and references therenn). .., , ,. t . . , . , . tU . , c 

C3 ™ 7- — , , , , , j. tidal disruption rate v, can be approximated to within an order of 

The capture radius of a black hole - the maximal distance .. , c ' , ■ / . , . , , , , , 

r ■ »;t,.. , •, ,r , , „ magnitude from an analysis of star dynamics near a black hole 

from the MBH where the tidal forces can overwhelm the stellar . ^ , , , ^ r . .. 

, r ■ , , ■ ■ via the Fokker-Planck equation. For the parameters of the GC it 

self -gravity and tear the star apart - is given by . . in-4 -l / lU • c m m 

b J t- b j gives the rate v s ~ 10 years (see the review of Alexander, 

/ M \'^ 3 /M I R \ 200 5J), which is in a greement with more detailed calculations of 

R T * 1.4 X 10 13 6 bh — — cm, (1) bver & UlmeJ (ll999h who obtained v f ~ 5 x 10 5 years" 1 . Tidal 

\10 M e ) \M Q ) \Rq) disruption processes were perha ps already observe d in cosmo- 

where M. and R t are the star's mass and radius, M bh the mass lo 8 ical § alax y surve y s ( see ' e.g.| Donley et al lEooi. 
of the black hole, and M and R Q the solar mass and radius. As observations show, there is a supermassive black hole 

When a star comes within the capture radius of a black hole, ( S f A *) in the center of our Galaxy with a mass of (4.3 1 +0.06)x 

the tidal force produc ed by the black hole i s strong enough to 10 M ° dGillessen et al. [ |200J). In the close vicinity of the 

capture the star JleTiL [1988; Phin^3, [T^9h . If the mass of the Galactic black hole (less than 0.04 pc from it) about 35 low mass 

black hole is higher than ~ 10 8 M Q , the star directly falls into the stars (1-3 M B ) and about lOmassive stars (3-15 M ) are present 

black hole event horizon and not much interesting phenomena ( see l Alexander & Liviol |2004D, i.e., processes of star accretion 

can be observed. However, if the mass of the black hole is less ln the Galactic center seems to be quite possible. However, at- 

than ~10 8 Mg, the star is torn apart. Roughly half of the star is tem P ts t0 find a stron S X " ra y source in the Galactic Center (GC) 

captured by the black hole and the disrupted debris form in a were falled - ^ Chandra X-ray Observatory resolved only a 

weak X-ray point sour ce at the position of S gr A* with a flux 

* cmko@astro.ncu.edu.tw L x ~ 10 33 erg s" 1 (see Baganoff et all |200"3|) . though moderate 
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X-ray flares were obser ved by Chandr a dPorauet et al.. 2008). 
As it w as mentioned by iKovama et al.l (1 19961) and iMuno et ail 
(2004) this "X-ray quiet" Sgr A* is in sharp contrast to the high 
X-ray activity of the surrounding diffuse hot plasmas. One likely 
scenario is that the Galactic nucleus was brighter in the past, 
possi bly caused by a surge ac cretion onto the massive back hole. 

ICheng et alj d2006l [2007) suggested that this scenario may 
explain the origin of the 511 keV annihilation flux from the GC 
region, if up to 10% of captured stellar matter is ejected in the 
form of a jet of relativistic protons. In this case the origin of the 
511 keV line emission from the GC region is supposed to be 
due to annihilation of secondary positrons generated by p - p 
collisions. 

An attempt to find independent evidences for active pro- 
cesses in the GC in th e form of a flux o f fast charged parti- 
cles was undertaken by Dogi el et~aT1 ((2008), who concluded that 
these primary relativistic protons (if were generated) penetrat- 
ing into dense molecular clouds produced there a flux of nu- 
clear de-excitation gamma-ray lines coming from the GC region. 
However, the origin of relativistic protons is still rather specula- 
tive since we do not have direct evidence s in favour of their pro- 
duction near black h oles (see, however, lAbraham et all [2007; 
IIstomin& Sol 12009) . Besides, the gamma-ray line flux would 
be strongly time variable in this case and, therefore, the estimates 
of its value at present are highly model-dependent. 

Below we examine a different model of gamma-ray line 
emission, assuming nuclear interactions of subrelativistic pro- 
tons generated by processes of star disruption at the GC. 

2. Flux of subrelativistic protons generated by star 
disruption 

An alternative (and to our view more reliable) mechanism of 
gamma-ray line production follows from the analysis of pro- 
cesses of star disruption near black holes. The energy budget of 
a tidal disruption event follows from analysis of star matter dy- 
namics. Once passing the pericenter, the star is tidally disrupted 
into a very long and dilute gas stream. The outcome of tidal dis- 
ruption is that some energy is extracted out of the orbit to unbind 
the star and accelerate the debris. Initially about 50% of the stel- 
lar mass becomes tightly bound to the black hole , while the 
remainder 50% o f the stellar mass is forcefully ejected (see, e.g. 
lAvalet aUl2000l) . The kinetic energy carried by the ejected de- 
bris is a function of the penetration parameter b^ 1 and can signif- 
icantly exceed that released by a norm al supernova (~ 10 51 erg) 
if the orbit is highly penetrating (see Alexander, 2005), 
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Thus, the mean kinetic energy per escaping nucleon is given by 
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where tjM„ is the mass of escaping material. For the black-hole 
mass M bh = 4.31 x 10 6 M the energy of escaping particles is 
£esc ~ 68(?7/0.5r I (fr/0.1)~ 2 MeVnucleon -1 when a one-solar 
mass star is captured. The parameters rj and b are not well con- 
strained by theory. It is clear, however, that b must be less than 
unity for a star to be disrupted in the gravitational potential of the 
black hole. We take the energy distribution of the erupted nuclei 
as a simple Gaussian distribution 
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where N is total amount of particles ejected by one stellar cap- 
ture. Results of calculations for different values of cr do not differ 
much from each other as long as cr < E esc (see below). 

For a single capture of a one-solar mass star the total number 
of unbounded particles is N i 10 57 and the total kinetic energy in 
these particles is W ~ 10 53 erg, if their energy is about 68 MeV 
(see above). For the star capture frequency v s — 10" 4 year -1 and 
E esc about several ten MeV (see Dogi el et al.l (12009b)) it gives a 
power input W ~ 3 x 10 41 erg s _1 , which is a few times more than 
the total power contained in Galactic cosmic rays. The average 
energy output depends, of course, on the capture frequency. 

A capture of a massive star should eject more fast particles 
into the GC region than a capture of a low mass star. However, 
the energy of the ejected particles is expected to be similar. This 
can be seen from Eq. (01, recalling that in first approximation 
the radius of main-sequence stars varies as R* oc M° 8 . The fre- 
quency of mas sive star capture is lo wer than that for low mass 
stars (see, e.g. ISver & Ulmerl 1 19991) . Therefore, we expect the 
capture of a massive star to have relatively little effect in com- 
parison with the cumulative effect from low-mass star captures. 
The latter is calculated below assuming M„ = 1 M for all cap- 
tured stars. 



3. Proton propagation in the Galactic center region 

Details of the proton spatial distribution in the GC region, as 
well as the mechanism of proton propagation are inessential for 
calculations of the total gamma-ray line flux. We use, neverthe- 
less, kinetic equations in which processes of propagation are in- 
cluded, and then integrate the obtained solution over the volume 
of emission. The effect of propagation reduces in this case to a 
proper estimate of the proton flux leaving the emission region. 

Propagation of cosmic rays in the Galaxy is described as a 
diffusion phenomenon, with a phenomenological diffusion coef- 
ficient whose value and energy dependence is derived from ob- 
servational data, e.g. from the chemical composition of cosmic 
rays measured near Earth. This method leads to a diffusion co- 
efficie nt of about 10 27 cm 2 s _1 (for details see lBerezinskii et all 
1990). Processes of particle propagation in the GC region are 
questionable since we do not know much about the physical con- 
ditions of the ambient medium. 

From the general theory of cosmic ray origin it is clear that 
particle propagation in the interstellar medium is described as 
diffusion due to scattering on magnetic fluctuations. The effec- 
tive coefficient of diffusion is D ~ v 2 /v(E,r), where v is the 
particle velocity and v(E, r) is the scattering frequency of par- 
ticles, which is a function of particle e nergy E and coor dinates 
r. Rather simple estimates provided by [Jean et al.l (120061) for the 
Galactic center region in the frame of the quasi-linear theory 
gives D ~ 10 27 cm 2 s' 1 for MeV pa rticles. A similar estimate 
was obtained by LaRosa et al.l d200ll) from the analysis of radio 
emission from the central region. These authors concluded that 
electrons of 0. 14 GeV diffuse about 13 pc in 4x 10 4 years, which 
gives again D ~ 10 27 cm 2 s _1 . 

Estimates of the diffusion co efficient for subrelati vistic pro- 
tons in the GC were derived by iDogiel et al.l (|2009c) from the 
observed hard X-ray emission, which was assumed to be due 
to inverse bremsstrahlung. They estimated D to be in the range 
10 26 - 10 27 cm 2 s _1 . Below we provide our calculations for the 
diffusion coefficient D = 10 27 cm 2 s _1 , though, in principle, its 
value may be a function of particle energy (as well as of spatial 
coordinates and even time). We suppose, however, that this sim- 
plification is acceptable because the energy range of relevance 
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for gamma-ray line production is quite narrow, less than one or- 
der of magnitude. 

Another important parameter of the model is the mean 
density of the medium into which the fast nuclei propagate. 
Observations show that the interstellar medium near the GC is 
highly nonuniform. Thus, within the central 2 pc of the shell of 
the supernova remnant Sgr A East, there is a hot p lasma of tem- 
peratu re 2 keV and density as high as 10 3 cirT 3 dMaeda et all 
|2002|) . This local density fluctuation is insignificant for our 
treatments since energetic particles capable of producing de- 
excitation gamma-ray lines lose little energy passing through 
this local fluctuation by diffusion. 

A significant part of the gas in the GC region is in the form 
of molecular hydrogen, ~9Q% of which is contained in dense 
and massive molecular clouds with a volume filling factor of 
only a few per cent. The total mass of molecular hydrogen in the 
Nuclear Bulge, which is estimated to be 2 x 10 7 M Q , significantly 
exceeds the mass con tained in the intercloud hot plasma (see 
lLaunhardt etaTll2002h . 

However, subrelativistic charged particles may not be 
able to penetrate freely insid e molecular clouds (see, e.g., 
Stalling & Str ong. A. Wlll976l) . Theoretical investigations per- 
formed by Dogiel et al. (1987, 2005) showed that turbulent mo- 
tions of neutral gas could excite small scale electromagnetic fluc- 
tuations that could prevent charged particles to penetrate deeply 
into clouds. Thus, subrelativistic cosmic-ray particles are ex- 
pected to fill only a small portion of molecular clou ds. This con- 
clusion is supported by calculations performed by Dogie lTtaT] 
(2009a), who showed that the flux of bremsstrahlung emission 
produced by subrelativistic protons inside molecular clouds is 
lower than that produced in the intercloud hot plasma. If this is 
true, we expect that most of the gamma-ray line emission should 
also be generated in the intercloud hot plasma. 

Recent observations with the satellites ASCA, Chandra and 
Suzaku showed that the 1-2° radius central region is filled with 
a hot gas of temperature 6-10 keV. The density of plasma de- 
rived from these observations rang e s between n ^ 0. 1 - . 4 cm" 3 
dKovama et all 119961; iMuno et all 12004 iKovama et all 120071) . 
In the calculations we use the average value n = 0.2 cm -3 . For 
these parameters the aver age time of Coulom b losses for 10-100 
MeV protons is (see, e.g. lHavakawa L H~969 ) 



E p m e w p 
6mie 4 In A 



3 ■ 10 6 - 3 ■ 10 7 yr, 



(5) 



where E p and \ p are the proton energy and velocity, and In A is 
the Coulomb logarithm. 



4. Spatial and energy distributions of subrelativistic 
protons in the Galactic center region 

The time-dependent spectrum of subrelativistic protons erupted 
from the central region can be calculated from the equation 
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where dE/dt = b(E) is the rate of proton energy losses. 
Subrelativistic protons lose their energy by Coulomb collisions. 
The rate of losses for these protons can be presented in the form 
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Fig. 1. Energy spectra of primary protons at r — pc and t = 
10 3 years {thick dashed line), r = pc and t = 10 4 years {thick 
solid line), r - 150 pc and t = 10 4 years (thin solid line), and r = 
300 pc, t — 10 4 years (thin dashed line). For r > 10 pc the proton 
spectra are almost stationary. These spectra were calculated for 
£ esc = 100 MeV and cr = 0.05£ esc - 



The injection of protons by processes of star capture can be de- 
scribed by 



Q(E,r,t)=Y J Qk(E)S(t-t k )S(r), 



(8) 



k=0 



where t k — kx T is the injection time and the functions Qk(E) are 
given by Eq. ©. The mean time of star capture by the massive 
black hole in the Galaxy is taken to be T - 10 4 years. 

We can derive the Green function of Eq . (0 for the injection 
time t k (see IS vrovatskiil [19591; iGrattonl 11972b and then sum over 
injections. The equation for the Green function is 

— - - V (DVG*) + — (b(E)G k ) = 5(E - E )6(t - t k )6(r) . (9) 
at oE 
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the Green function can be written as 
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In the nonrelativistic case, i.e. for energy losses in the form 
of Eq. (0, the solution ( fT2l can be simplified to 
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Fig. 2. Densities in protons of energies 10, 50, 80 and 100 MeV, 
as a function of radial distance from the GC. These spatial distri- 
butions were calculated for E esc = 100 MeV, <x = Q.Q5E esc , and 
n - 0.2 cirT 3 . 



As one can see this solution is characterized by two different 
spatial scales. The first one is R\ ~ V4DT ~ 10 pc. Inside the 
sphere r < R\ the solution is time dependent and the density of 
protons is strongly fluctuating with time. On the other hand, for 
r > R\ the solution is quasi-stationary because the time of diffu- 
sion over these distances is larger than the mean capture time T. 
The radius of the sphere filled with subrelativistic protons (the 
scale of exponential decrease) is about Ri ~ -\f4Dri ~ 200 pc 
for £gsc = 100 MeV. One can see from Eq. (0 that the mean 
free path of protons is proportional to y/l/n. Examples of cal- 
culated energy and spatial distributions of subrelativistic protons 
are shown in Figs.[T]and|2] 

Since the Galactic black hole was inactive for a long time, 
we use to calculate the gamma-ray line emission the proton spec- 
trum at time t — T after the last star capture. This will give us 
a lower limit on the gamma-ray line flux. However, the density 
of subrelativistic protons in the GC is expected to be almost in- 
dependent of time, except for a very compact region around the 
massive black hole. The proton spectrum should reach its satu- 
ration level after > r,/r capture events. In our calculations we 
summarized the cumulative effect of 2000 captures. 

In order to predict the total gamma-ray line flux from the GC 
we need to estimate the total (integrated over volume) spectrum 
of subrelativistic protons 



F(E,t) = 4ttJ f(r,E,t)r 2 dr. (15) 

o 

We note that for simple estimates one can derive this spec- 
trum from the assumption of a stationary injection of protons 
with the rate Q = N/T. Then the approximate solution has the 
form 

F{E) = \\ +Erf ( Eesc ~ E )] (16) 

2a [ \ V2cr /. 

where Erf(x) is the error function. 



5. Nuclear interaction gamma-ray line emission 
from the Galactic center region 

Unlike relativistic particles, subrelativistic nuclei do not pro- 
duce effectively continuous radiation (except, may be, inverse 
bremsstrahlung, which is discussed in Dogiel et al. 2009abc). 
However, collisions of subrelativistic nuclei with ambient mat- 
ter can lead to nuclear excitation and result in emission of de- 
excitation gamma-ray lines. These lines may be a good tracer 
for subrelativistic cosmic rays, because the line brightness can 
give us information about the amount of subrelativistic particles. 

Thus, prominent nuclear de-excitation lines are observed in 
emission spectra of solar flares, which allows one to derive valu- 
able information on solar ambient abundances, density and tem- 
perature, as well as on accelerated particle composi tion, spec- 
tra and tra nsport in the solar atmosphere (see e.g.. ISmith et al.l 
(2003]) and lKiener et al.l d2006l) for rece nt solar observa tions with 
RHESSI and INTEGRAL, respectivelv). |Pogiell d200ll) predicted 
that galaxy clusters could emit a detect a ble flu x of de-excitation 
gamma-ray lines. Latter, llvudin et al.l d2004l) found traces of 
gamma-ray line emission towards the Coma and Virgo clusters, 
though this detection has not been confirmed yet. 

Detection of the gamma-ray line emission produced by 
cosmic-ray interactions in the interstellar medium would provide 
insightful information on low-energy cosmic rays and give a sig- 
nificant advance in the development of the theory of cosmic -ray 
origin. Theoretical estimates for the lin e emission in the Galaxy 
were provided by Ram atv et al.l (1 19791) . who calculated spectra 
of the nuclear de-excitation line emission for different assumed 
spectra of subrelativistic ions. However, attempts to measure a 
prominent flux of nuclear interaction lines were failed up to now. 

Galactic center observations in the MeV range were per- 
formed by the COMPTEL group but with a rather poor angular 
resolution dBloemen et all Il997l) . They detected a marginal ex- 
cess of emission in the range 3-7 MeV. The total flux from the 
central region 60°x20° was estimated to be 10~ 4 cm~ 2 s _1 rad -1 . 
This excess was interpreted as emission in nuclear de-excitation 
lines generated by cosmic-ray interaction in the GC region. 
However, neither the COMPTEL group nor subsequent obser- 
vatio ns with INTEGRAL have con firmed this result. In particu- 
lar, |^egarden&^^tan^3 (E006) have carried out an extensive 
search for gamma-ray lines in the first year of public data from 
the INTEGRAL spectrometer SPI. They found no evidence for 
any previously unknown lines in the energy range 20-8000 keV. 

We suppose that processes of star accretion by the mas- 
sive black hole at the GC can generate a significant number of 
subrelativistic protons and heavier nuclei. These energetic par- 
ticles could in turn produce a significant flux of nuclear de- 
excitation gamma-ray lines in a relatively compact region of the 
GC. Unlike the ga mma-ray emis s ion th at might be produced by 
relativistic protons dDogiel et al.L [2008). the emission produced 
by subrelativistic protons is expected to be almost stationary and 
should therefore be less dependent on model parameters. 

Though according to the present model the origin of subrel- 
ativistic cosmic rays near the GC is different than in other part 
of the galactic disk, a positive detection of de-excitation gamma- 
ray lines would provide important information about star accre- 
tion processes at the GC. 

We calculate the gamma-ray flux emitted at time t — T after 
a star capture in a given nuclear de-excitation line from 

F rW = j-^y— FtiEuOviiEdo-tjiEddE,, (17) 
And 2 *rji n H Jo 
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Fig. 3. Radial profiles of enclosed fluxes in (a) the 4.44 MeV 
line and (b) the ~6.2 MeV line complex, as a function of angular 
radius from the GC, for the five values of £ e sc indicated in the 
figure. The calculations assume cr = Q.Q5E esc (Eq. |4ji and «h = 
0.2 cirT 3 . 



where i and j range over the accelerated and ambient particle 
species that contribute to the production of the gamma-ray line 
considered, «h = 0.2 cirT 3 is the adopted mean density of H 
atoms in the interaction region, rij is the density of the ambient 
constituent j, Fi{Ei, t) is the energy spectrum at time t of the fast 
particles in the interaction region (Eq. 15), Vj is the velocity of 
these particles, cry is the cross section for the reaction of inter- 
est between species i and /, and d - 8 kpc is the distance to 
the Galactic center dGronewegen et all 120081) . To take into ac- 
count the measured en hancement of the m etal abundances in the 
Galactic center region (ICunha et al.Ll2007l) . we assume the abun- 
dances of ambie nt carbon and heavier elements to be twice solar 
dLoddersl 12003). This is consistent with the recent estimate from 
Suzaku observations of[kovama et al. (2009), who found the Fe 
abundance to be between 1.5 and 2.3 times solar. In the case of a 
low-mass star capture, we expect the line emission to be mainly 
due to proton and cc-particle reactions with ambient heavy nu- 
clei. We took into account the a particles assuming, for simplic- 
ity, F a (E, f) - X a F p (E, t), where the energy E is expressed in 
MeV nucleon -1 and X a — 0.1 is the accelerated a/p abundance 
ratio. 

Noteworthy, the total steady-state emission obtained from 
Eq. ( TP7T i is independent of the assumed mean density in the 
interaction region, as can be seen by introducing in this equa- 
tion the particle distribution Fj(E, t), which is inversely propor- 
tional to «h (see, e.g., Eq. [Tol l. This is true, however, as long as 
«h ~ 100 crrT 3 , because for larger density values the Coulomb 
loss time t,- < T and the particle distribution is not anymore sta- 
tionary. In any case, even if the predicted emission were time de- 
pendent, the total fluxes given below can be interpreted as mean 
values over the period of star capture T. 

Figure|3]shows calculated radial profiles of the predicted line 
emission at 4.44 and ~6.2 MeV. These two lines are mainly pro- 
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Fig. 4. Predicted fluxes from the GC region of angular radius 
9 = 5° in (a) the 4.44 MeV line and (b) the -6.2 MeV line 
complex, as a function of i? e sc, for cr = 0.05 £ e sc (solid lines) and 
cr = E esc (dashed lines). 



duced by proton and o--particle reactions with ambient 12 C and 
16 (see below). Contrary to the total emission from the GC re- 
gion (i.e. integrated over space), the radial profiles depend on the 
assumed value of nu' the larger the ambient medium density the 
lower the spatial extend of the gamma-ray line emission. We see 
in Fig.|3]that for «h = 0.2 cirT 3 and E esc < 100 MeV nucleon -1 , 
the gamma-ray lines are emitted from a region of maximum 
5° angular radius. Noteworthy, such an emission would appear 
as a small scale diffuse emission for a gamma-ray instrument 
like SPI, whose angular reso lution is ~3° FWHM (see, e.g., 
iTeegarden & Watanabd 120061 

Figure [4] shows calculated fluxes in the 4.44 and 
~6.2 MeV lines as a function of E esc ranging from 10 
to 100 MeV nucleon 1 , which corresponds to 0.68 < 
(?7/0.5)(£/0.1) 2 < 6.8 (see Eq. |3). We see that the predicted 
fluxes increase with increasing escape energy up to values « 
10~ 5 photons cirT 2 s -1 for E &sc * 100 MeV nucleon -1 . We also 
see that the fluxes are not strongly dependent of the assumed 
wi dth of the adopted Ga ussian distribution (cr in Eq. 4). 

DogieTetaTJ (|2009b) showed t hat the hard X-rays observed 
with Suzaku from the GC region dYuasa et"aTJ, 12008b could be 
produced by inverse bremsstrahlung from subrelativistic protons 
accelerated by star disruption near Sgr A*. In this scenario, the 
measured spectrum of this hard X-ray emission would imply 
that E esc « 100 MeV nucleon -1 . Figure shows a gamma-ray 
line spectrum calculated for E esc = 100 MeV nucleon" 1 an d 
cr = 0.05 £esc, using the code developed bv lRamatv et alJ dl979). 
The calculated total flux in gamma-ray lines below 8 MeV 
is 1.1 x 10 -4 photons cm -2 s -1 . In the 3-7 MeV range, it is 
4.4 x 10 -5 photons cm -2 s -1 . 

We see in Fig. [5] relatively strong lines of various widths, 
which are superimposed on a continuum-like emission due to 
a large number of unresolved gamma-ray lines, mostly aris- 
ing from cascade transitions in high-lying levels of heavy nu- 



6 



Dogiel et al.: Nuclear interaction gamma-ray line emission 




1 
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Fig. 5. Spectrum of the predicted gamma-ray line emission from 
the GC region of angular radius 9 - 5°, assuming E esc 
1 00 MeV nucleon 1 and cr = 0.05£esc- The arrows point to the 
lines which are mentioned in the text. 



clei. Among the strongest lines are those produced in excita- 
tions of low-lying nuclear levels of abundant ambient ions, e.g. 
at 4.44 MeV from 12 C* and 0.847 MeV from 56 Fe*. Other im- 
portant lines are those arising from transitions in the spallation 
products of abundant ambient nuclei, e.g. at 0.718 MeV from 
10 B*, which is produced by spallation of ambient I2 C and 16 0, at 
0.931 MeV from 55 Fe*, which is mainly produced by the reac- 
tion 56 Fe(j?,/?«) 55 Fe*, and at -5.2 MeV from 15 N* and 15 Q * re- 
sulting from spallation of ie O (see iKozlovskv et ail, 12002). The 
broad line feature at -6.2 MeV is due to the merging of the 
6.129 MeV line from de-excitations in ie O* with two lines at 
6.175 MeV from ls O* and 6.322 MeV from I5 N* produced in 
spallation reactions. The characteristics of the main gamma-ray 
lines are summarized in Table [TJ 

The lines produced in a gaseous ambient medium are broad- 
ened by the recoil velocity of the excited nucleus. The line 
width generally decreases with increasing mass of the target 
atoms, because the recoil velocity of heavier nuclei is smaller. 
Thus, in Fig. [5] the full width at half maximum (FWHM) of 
the 0.847 MeV line produced by proton and or-particle inelas- 
tic scattering of 56 Fe is 6 keV (0.7% of the transition energy), 
whereas the FWHM of the 4.44 MeV from 12 C* is 160 keV 
(AEy/Ey = 3.6%). 

However, some gamma-ray lines produced in interstellar 
dust grains can be very narrow, because some of the excited 
nuclei can stop in solid materi als before emitting gamma-rays 
dLingenfelter & Ramatvl 1 19771) . It requires that the mean life- 
time of the excited nuclear level or of its radioactive nuclear par- 
ent is longer than the slowing down time of the excited nucleus 
in the grain material. The most promising of such lines is that at 
6.129 MeV from ie O*. Thus, future measurements of the profile 
of the -6.2 MeV line could allow us to probe the presence of 
dust grains in the Galactic center region and perhaps t o measure 
the grain size distribution (Tatischeff & Kienerl 12004). 

Very narrow lines can also result from the decay of 
spallation-produced long-lived radionuclei, which can come es- 
sentially to rest in the ambient medium before decaying to ex- 
cited states of their daughter nuclei and emitting gamma-rays. 



Table 1. Main gamma-ray lines expected from the GC region 



Line energy Nuclear FWHM" Flux" 

(MeV) Transition (keV) cirT 2 s _1 



0.718 


10 R « . 
°0.718 l 


>.S.' 


30 


3.4x10-" 


0.847 


rc 0.847 


g.S. 


6 


1.4xl0~ 6 




27 A1* 
rtl 0.844 


g.S. 






1.809 


ivl 5l 8(19 ^ 


g.S. 


<2.7 C 


3.5xl0~ 7 




( 26 A1 decay) 






4.44 


12 C* 


g.S. 


160 


1.2xl0~ 5 




Ilr>« . 
D 4.445 


g.S. 






-5.2 


14 N * _^ 
^5. 106 


g.S. 


280 


5.4xl0" 6 




15 o* -> 

w 5.181 


g.S. 








15 o* 

w 5.241 ^ 


g.S. 








^5.270 


g.S. 








y 5.299 


g.S. 






-6.2 


16 o* 

w 6.130 


g.S. 


180" 


1.0xl0~ 5 




15 o* 


g.S. 








^6.324 


g.S. 







" For £ es c = 100 MeV nucleon" 1 and a-=0.05£ esc (see Eq. 4). The 
fluxes are for the region of angular radius 6 = 5° (see Fig. [3] for 
fluxes at 4.44 and -6.2 MeV calculated for lower values of ff). 

b Ground state. 

c Assuming thermal broadening with T < 10 keV. 

d A very narrow line component at 6. 129 MeV could arise from 
excitation of ls O nuclei contained in micrometer-sized grains (see 
text). 

For «h ~ 0.2 crrT 3 , the slowing-down time of recoiling heavy 
nuclei excited by fast proton and c-particle collisions is of the 
order of 10 4 years. Thus, only radionuclei with mean lifetimes 
longer than 10 4 yr can produce very narrow lines in the GC envi- 
ronment. The most intense of these lines is expected to be that at 
1 .809 MeV from the decay of 26 Al (70 1/2 = 7.2 x 10 5 years). This 
radioisotope can be synthesized by the reactions 26 Mg(p,n) 26 Al 
and 28 Si(p,Jt) 26 Al (iRamatv et all 1 19791) . However, the calculated 
flux in the 1.809 MeV line is only 3.5xl0" 7 photons crrT 2 s -1 
(Table[TJ. In comparison, the flux observed with INTEGRAL/SPl 
from the inner Galaxy (-30° < I < 30°; - 10° < b < 10°) is 
(3.3 + 0.4) x 10~ 4 crrT 2 s _1 rad~'. It is produced by a present-day 
equilibrium mass of 2.8 + 0.8 M Q of 26 Al s ynthesized by massive 
stars in the Milky Way (iDiehl et al .1120061) . 

The gamma-ray lines above 4 MeV are better candidates for 
a future detection. Our predicted fluxes, however, are below cur- 
rently available sensitivity limits of INTEGRAL. After one year 
of observation and -3 Ms of exposure of the GC region, the 
SPI sensitivity for detection of a narrow line at -5 MeV from 
a small-scale diffu se region at the center of the Galaxy was 
3.2xl0~ 5 cnr 2 s _1 iTeegarden & WatanabeLl2006l see Table 5). 
But for a broad line of -200 keV FWHM (see Table Q} the sen- 
sitivity limit increases by a factor of -8. It is thus unlikely that 
the predicted gamma-ray line emission will be detected with 
INTEGRAL. 



6. Discussion and conclusions 

In spite of many uncertainties of the model parameters we should 
state that their estimates cannot vary more than an order of mag- 
nitude. Thus, the injection energy should be close to 100 MeV, 
the average energy output is about 10 41 erg s _1 and the spa- 
tial diffusion coefficient is about 10 27 cm 2 s _1 in order to repro- 
duce the observed spatial distribution of hard X-rays in the GC 
(Dogiel et al. 2009a,b,c). 
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We showed that accretion of stars onto the central supermas- 
sive black hole periodically eject in the GC region an intense 
flux of subrelativistic protons and nuclei capable of producing 
a significant emission of nuclear de-excitation gamma-ray lines. 
The production rate of these energetic particles is of the order 
of 10 45 s _1 . Because the time period of star capture by the black 
hole is expected to be shorter than the lifetime of the ejected 
protons against Coulomb losses, the density of these particles in 
the GC region is expected to be almost stationary. The radius 
of the volume filled with these energetic particles is found to be 
i 700 pc or t 5° if observed from Earth. 

Based on the spectrum of the diffuse h ard X-ray emissio n 
observed with Suzaku from the GC region dYuasa etall [2008), 
we estimated the energy of the nuclei accelerated in the vicinity 
of the black hole to be E esc » 100 MeV nucleon -1 . Further as- 
suming that the mean metallicity in the GC regi on is two times 
highe r than in the solar neighborhood (see, e.g., iKovama et al.L 
2009), we calculated the total gamma-ray line flux below 8 MeV 
to be 1 . 1 x 10~ 4 photons crrT 2 s _1 . The most promising lines for 
detection are those at 4.44 and ~6.2 MeV, with a predicted flux 
in each line of »T0~ 5 photons cirT 2 s _1 . These lines should be 
broad, AE 7 /E y of 3^4%, which unfortunately renders their de- 
tection with the INTEGRAL spectrometer unlikely. 

But future ga mma-ray mi s sions like the Nuclear 
Compton Telescope dChang et all 120071) . the GRIPS project 
(Greiner J., et al. , 2008), and the Advanced Compton Telescope 
(Bog gs et all |2006j) may be able to test our predictions. In 
particular, the GRIPS mission proposed for ESA's "Cosmic 
Vision" program could achieve after 5 years in orbit more than 
an order of magnitude sensitivity improvement over COMPTEL 
(in 9 years), which would allow a clear detection of the predicted 
gamma-ray line emission at 4.44 and ~6.2 MeV from the GC 
region. The Advanced Compton Telescope project proposed 
as a future NASA mission aims at even better sensitivity, near 
10~ 6 photons cirT 2 s _I for 3% broad lines. A future detection 
of the predicted gamma-ray lines with such an instrument 
would provide unique information on the high-energy processes 
induced by the the central black hole, as well as on the physical 
conditions of the emitting region. 

Finally, we note that nuclear interactions of subrelativistic 
ions with ambient material can also synthesize f3 + radioisotopes, 
whose decay can inject positrons into t he GC region. From th e 
radioisotope production yields given bv lKozlovskv et all d 19871) . 
we estimate that for E esc a; 100 MeV nucleon 1 the number of 
positrons produced by this mechanism is i 5 times the number 
of gamma-rays emitted in the 4.44 and ~6.2 MeV lines. With a 
predicted gamma-ray production rate of 1 .7 x 10 41 photons s _1 in 
the sum of these two lines (see Table 1), it gives a positron pro- 
ductionrate i 8.4xl0 41 y6 + s . This limit is more than an order 
of magnitude lower th an the positron annihilation ra te measured 
with INTEGRAL/ S PI dWeidenspointner et aUl2008h . 
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